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Abstract 

Recently the AGASA Collaboration presented data suggesting a significant clus¬ 
tering of ultra-high energy cosmic rays coming from the outer Galaxy region. In 
this paper we calculate expected cosmic ray arrival distributions for several simple, 
limiting source location scenarios and investigate the possibility of clustering and 
correlation effects. The role of the Galactic magnetic field is discussed in detail. 


1 Introduction 

The observations of cosmic rays of energy above 10^® eV reveal at least two features that 
we have not yet understood and that appear to contradict each other. The cosmic ray 
energy spectrum does not seem to be cut off because of photoproduction interactions 
on the microwave background 0, H and extends above 10^° eV [§]. The existence of 
these particles suggests that the sources of ultra-high energy cosmic rays (UHECRs) 
are cosmologically nearby - within 20 Mpc or so @]. On the other hand, the UHECR 
arrival direction distribution, as far as we know it, seems isotropic on a large scale with a 
statistically improbable small scale clustering p (for a review on the subject see [^). 

This small-scale clustering was hrst observed in the data of the AGASA group |^. The 
current AGASA data set of 59 events of energy above 4 x 10^® eV contains hve doublets 
and one triplet for a maximum separation angle of 2.5° P]. The centroids of the triplet 
and one of the doublets are less than 1 degree off the supergalactic plane (SGP). When 
combined with data from the Haverah Park, Yakutsk and Volcano Ranch experiments 
the number of doublets (within 3°) increases to eight and the number of triplets to two 
ip. The chance probability of observing these multiplets from a uniform distribution of 
sources is less than 1.5 % and less than 1 % when a restricted region within ±10° off the 
SGP is considered [Q. 
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Recently the AGASA group presented a self-correlation analysis of their data with 
E > 10^® eV [|]. Remarkable correlations are found in Galactic coordinates, supporting 
further the previously found indications of clustering of UHEGRs. In a straight-forward 
interpretation one would link the existence of such clusters directly to UHEGR point 
sources |P, |T^ but non-trivial effects such as possible clustering of sources or large-scale 
magnetic helds could also contribute to such a correlation |^, |^, |^, . 

The indications found in data raise the interesting question of whether both an isotropic 
distribution and a clustering of the arrival directions can be explained consistently, as¬ 
suming they are not a statistical fluctuation [^, 

Potential models of UHEGR origin, assuming they are charged particles, are con¬ 
strained by their ability to reproduce the measured energy spectrum and the approxi¬ 
mately isotropic arrival distribution observed in the data [PU, 




At the same time, their 
predictions should also be consistent with the small scale clustering. In general, several 
key ingredients enter the models: (i) the distribution (locations) of the sources, (ii) their 
nature, i.e. whether they are sources emitting cosmic rays continuously for a long period 
of time or they are bursting sources, (iii) the energy spectrum and total flux at injection, 
and (iv) the propagation of the cosmic rays including energy loss processes and deflection 
due to magnetic helds. Glustering on small scales and overall isotropy seems to favour 
models based on several nearby sources (for example, in the Galactic halo). Models of this 
type include the acceleration of cosmic rays at magnetars m , UHE neutrino interactions 
on the relic neutrino background (Z-bursts) [P3, 1^, and the decay and annihilation of 


superheavy Relics accumulated in the halo of the Galaxy |^, |2^, |2^, ^ 


In this paper we study general aspects related to cosmic ray propagation and the source 
distribution. The aim is to explore what a future conhrmation of the correlations could 
reveal about the nature of the sources of UHEGRs, their distribution and the strength of 
the magnetic held in our Galaxy. To keep our results as general as possible we consider 
three simplihed, limiting source location scenarios: 

(i) Uniform distribution of sources, which is relevant to UHEGR of extragalactic origin 
with sources of isotropic and homogeneous distribution. 

(ii) sources distributed uniformly within ±20° oh the supergalactic plane (SGP), and 

(iii) single and multiple point sources. 

For each of these source scenarios we calculate arrival distributions, one and two-dimensional 
correlation functions and, where applicable, arrival time delays. 

It should be emphasized that we do not attempt to reproduce AGASA data with 
these source models. The interest of this work is merely the investigation of characteristic 
features which are fairly model-independent. For a more realistic analysis one needs a 
better knowledge of the Galactic magnetic held, as well as full access to the experimental 
data, including acceptance corrections. 

The outline of the paper is as follows. In Sec. we discuss the Galactic magnetic held. 
The simulation method for cosmic ray arrival distributions is explained in Sec. ^ and the 
diherent source location scenarios are presented in Secs. ^ through |^. Sec. |] applies the 
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results of the previous sections to AGASA multiplets and Sec. concludes the paper with 
a summary of our main results. 


2 Galactic magnetic field 


Unfortunately the Galactic magnetic held (GMF) is not very well known. The position 
of the Solar system makes it difficult to measure its global structure. The most extensive 
and reliable information on the GMF comes from measurements of polarized synchrotron 
radiation and Faraday rotation of the radiation emitted from pulsars and extragalactic 
sources. The structure of the GMF as well as its local value at the Solar system are 
still uncertain, mainly due to the limited number of rotation measures available and 
the intrinsic difficulties in distinguishing local small-scale features from a large-scale one. 
Faraday rotation measurements indicate that the GMF in the disk of the Galaxy has a 
spiral structure with held reversals at the optical Galactic arms (for a review see [^). 
Recent work favours a bisymmetric spiral held (BSS) structure but an axisymmetric (ASS) 
held is not excluded [^, Q. 

In the following we use a BSS held model for the regular magnetic held in the disk of 
the Galaxy. The Solar System is at a distance rj| = = 8.5 kpc from the center of the 

Galaxy in the Galactic plane. The local regular magnetic held in the vicinity of the Solar 
System is assumed to be ~ 1.5 fiG in the direction I = 90° -|- p where the pitch angle is 
p = —10° @]. New measurements discuss larger total held strengths of up to 6 pG 
Therefore our assumptions should be considered as rather conservative. To illustrate the 
ehect of a stronger magnetic held we will also perform a simulation with a local regular 
held strength of 3 pG. 

The strength of the spiral held at a point in the Galactic plane, described by the polar 
coordinates (r||,0), is given by: 





cos ( (j) 


r\\ 
ro 


( 1 ) 


where Bq = 4.4 pG, ro = 10.55 kpc and /3 = 1/ tanp = —5.67. The held decreases with 
Galactocentric distance as l/r|| and it is zero for r|| > 20 kpc. In the region around the 
Galactic center (r|| <4 kpc) the held is highly uncertain. For simplicity we assume it is 
constant and equal to its value at r\\ = 4 kpc. Fig. ^ shows the held in the Galactic plane. 
For clarity we have not plotted the region r|| < 3 kpc. 


Following the spiral held strengths above and below the Galactic plane are taken 
to decrease exponentially with two scale heights: 


\B{r\\,(i),z)\ = \B{r\\,(i))\ 


exp(—z) 

exp(-3/8) exp(-z/4) 


\z\ < 0.5 kpc 
\z\ > 0.5 kpc 


( 2 ) 


where the factor exp(—3/8) makes the held continuous in The BSS spiral held we use 
is of even parity, i.e. the held direction is preserved at disk crossing. 


3 







X (kpc) 


Figure 1: Magnetic field configuration in the Galactic plane. The vectors indicate the 
held direction and their length is proportional to its magnitude. 


Observations show that the held in the Galactic halo is much weaker than that in the 
disk. The currently discussed options favour either a quadrupole held or a simple dipole 
held [^, 1^, In this work we assume that the regular held corresponds to a AO dipole 
held as suggested in |^. The dipole held is toroidal and its strength decreases with 
Galactocentric distance as 1/r^. In spherical coordinates (r, 9, ip) the {x, y, z) components 
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of the halo held are given by: 


Bx = —3 jjiQ sm.9 cos9 

By = —?) jjiQ sm. 9 cos 9 sin (3) 

Bz = /iG (1 - 3sin2 6')/r^ 


where ~ 184.2 /iG kpc^ is the magnetic moment of the Galactic dipole. The dipole 
held is very strong in the central region of the Galaxy, but is only 0.3 /iG in the vicinity 
of the Solar system, directed toward the North Galactic Pole. 

There is a signihcant turbulent component, i?rarn of the Galactic magnetic held. Its 
held strength is difficult to measure and results found in literature are in the range of 
-Bran = 0.5... 2i?reg [0 • To simulate it we add to both regular components a random 
held with a strength of 50% of the local regular held strength and a coherence length 
of 100 pc. The possible time dependence of the turbulent magnetic held component is 
neglected. This could be an over-simplihcation if this component changes signihcantly 
over time scales of 10^ years. 


3 Simulation and analysis methods 

3.1 Simulation 

We simulate the propagation of protons of energy above 10^® eV in the Galactic magnetic 
held (GMF) by numerically integrating the equations of motion in a magnetic held. The 
results can be easily re-scaled in rigidity for heavier nuclei. 

Since we are not interested in extragalactic propagation it is sufficient to consider a 
number of typical Galactocentric source distances, namely rsrc = 15, 20 and 40 kpc. The 
40 kpc distance corresponds to a scenario of distant halo sources or extragalactic sources 
if the extragalactic magnetic held is weak, and the two others are chosen to study the 
distance dependence within the Galaxy. The energy loss of protons can be neglected for 
these distances. Therefore we use a diherential injection spectrum which is similar to the 
observed spectrum, dN/dE cx 

Protons are injected at the corresponding source locations and each trajectory is fol¬ 
lowed until: 

(i) it reaches a Galactocentric distance larger than 100 kpc, 

(ii) the total path traveled by the proton is larger than 3 Vgrc, or 

(hi) it intersects a spherical detector of radius rdet centered at the position of the Solar 
System. 

As will become clear later, the total length of trajectories reaching Earth is always much 
smaller than the limit imposed by (ii). 

The radius of the detector is rdet = 100 (rsrc/20 kpc) in pc. This guarantees the 
same geometric efficiency for different source distances and optimizes the efficiency of the 
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propagation, keeping the nnavoidable smearing in arrival angle smaller than 0.5°. 

In all calculations we consider isotropically emitting sources. However, to increase the 
efficiency of the simulations, we only inject protons within a cone of half opening angle 
of 30°, pointing toward the detector. We have checked that the number of protons that 
reach the detector and are injected at an angle larger than 30° from the line-of-sight is 
negligible for all energies and distances considered here. 



Figure 2: Examples of proton trajectories projected onto the xz plane, being orthogonal 
to the Galactic plane. The left (right) panel corresponds to a source located at 6 = 
—45°, I = 112.5° {b = 0°, I = 112.5°) in Galactic coordinates at a distance of Vsic = 20 
kpc. The tracks are labeled with their corresponding energies. The dot marks the position 
of the detector. Note the different scales of the horizontal and vertical axes. Only the 
regular GMF has been used in this simulation. 

To illustrate our simulation procedure. Fig. ^ shows the projection of several represen¬ 
tative proton trajectories onto the xz plane. Only the regular GMF has been used in this 
simulation. The coordinate system is dehned by +z pointing toward the North Galactic 
pole and the location of the Solar system being at = (8.5, 0, 0) in units of kpc. The 
trajectories shown in the right panel clearly reflect the structure of the BSS field with its 
field reversals. 

3.2 Analysis 

One of the signihcant deviations from isotropy observed by the AGASA detector (and 
supported by the other air shower data sets) is the non uniformity in the distribution 
of space angle between the UHEGR events. This is caused by the large fraction of the 
experimental statistics above eV that is in the form of doublets - 13 out of 59 

individual showers. 
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The number of multiplets for a given opening angle depends on the total number of 
cosmic rays observed within the field of view of the detector. With increasing number 
of detected cosmic rays the number of “artihcial” multiplets grows simply due to the 
decreasing mean space angle between the cosmic rays. This is a particular problem for 
numerical simulations which are not restricted in statistics. Therefore we will, following 
the similar approach of AG AS A calculate the signihcance of one- and two-dimensional 
self-correlations. 

For each of the scenarios we calculate the self-correlation in separation angle of the 
events in the simulated sample. The separation angle between pairs of events is obtained 
and the distribution is divided by the solid angle of a concentric ring centered in each 
angular bin. The statistical signihcance of any deviation from the expected isotropic back¬ 
ground is shown in the signihcance plots as (signal-background)/sigma-of-background. 
The background is given by the mean self-correlation function for many sets with the 
same number of cosmic rays but arrival directions sampled according to a locally isotropic 
hux. 

The role of the Galactic magnetic held as well as possible bias due to cosmic ray selec¬ 
tion criteria can be studied by two-dimensional self-correlations in Galactic coordinates 
I, b. The diherence in Galactic longitude (AZ) and Galactic latitude (A6) between each 
pair of cosmic rays is calculated and plotted in a 2 dimensional map. The signihcance plots 
are meant to represent graphically the deviation of a certain arrival direction distribution 
from a local isotropic distribution. The function plotted is: 

p(Ai, Ab) = 

- {UJAl,Abr 

The background (bkg) is calculated by sampling from an isotropic distribution many sets 
of events containing the same number of events as the simulated data set. A large number 
of these sets has to be generated since the background depends strongly on the number 
of events falling into the A/ x Ab window of the correlation plot. Typically we generate 
10,000 such background conhgurations. 

The experimental resolution is important. We assume the error of a measured arrival 
direction is Gaussian. The distribution in arrival directions is then given by 

9('. i) = E A ^ 

T I W W 



where Uj is the error of the cosmic ray i with arrival direction (/j, bi). For simplicity we 
have assumed the same uncertainty in Galactic longitude and latitude. 

The correlation function reads 


= E 


/dat(AZ, Ab) = I dl I db g{l, b)g{l - A/,6 - Ab) 
(A/ - A/p)2 {Ab 


p ‘^TTO'idj 


exp 


Abp)^ 




crt: 


( 6 ) 

(7) 
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The sum runs over all pairs p = {i,j) of cosmic rays with i ^ j and A/p = k — Ij, 
Abp = bi — bj. The new effective error is 


(^P = \/(^i +(^], ( 8 ) 

as expected from standard error propagation. The plots are made with an energy- 
independent, constant error of cTj = 1.47°. Note that these plots are by construction 
symmetric, i.e. every pair of particles enters the plots twice. 

The signihcance plots can be interpreted as positive or negative deviation (over- or 
underdensity) in units of the uncertainty given by the expected background fluctuations. 
Consequently, the signihcance of an excess seen in a plot with only a few cosmic rays 
falling in the A/ x Ab window is always small. The large background huctuations reduce 
the signihcance of any possible signal. 

Since the particle gyroradius, and thus the degree of dehection in the GMF depends 
on its energy, the physics of the dehection of cosmic rays in the Galactic magnetic held is 
better revealed by looking at energy-ordered pairs. For example, knowing the sign of the 
charge of the cosmic rays, the energy-ordered correlation function can be used to derive 
the mean magnetic held direction along the path of weakly dehected cosmic rays coming 
from the same source. In the following we apply the convention 


A/ 



li Ij 

/j /^ 


Ei > Ej 
Ei < Ej 


(9) 


and the corresponding dehnition for Ab. For cosmic rays of the same charge, the corre¬ 
sponding correlation function is not symmetric in Al and Ab if there exists a magnetic 
held component orthogonal to the particle trajectories (see below). 

To check our analysis methods and to demonstrate the importance of the newly intro¬ 
duced quantities (Eqs. (^) and (^)) we analyse the published list of AGASA events with 
energy above eV, which consists of 58 events containing 1 triplet and 5 doublets 

3^. 

Following the AGASA analysis [|| we apply the cuts 90° < / < 180° and —60° < b < 
60°. This cut in arrival direction will be referred to as the “outer Galaxy” cut from now 
on. From the AGASA data 25 events with E > 10^®'® eV pass the outer Galaxy cut. They 
contain two of the doublets and the triplet. 

The left panel of Fig. shows the ID angular correlation of the AGASA data above 
1019.6 gY Qoming from the outer Galaxy. The significance of the deviation from an uniform 
distribution is plotted in the right panel. The signihcance approaches 5 ct for a space angle 
of 2° and is consistent with no deviation for larger space angles. 

Both the symmetrized and the energy-ordered signihcance of the self-correlation p{Al, Ab) 
of the same event sample are given in Fig. |[ Note that the signihcance of these 2D self¬ 
correlations is respectively 4.2cr and 4.6cr. The energy-ordered plot also shows the direction 
of dehection as a function of the energy. The signihcance peaks at Al ^ Ab ^ —1°. 







Figure 3: ID correlation function for published AGASA data with E > 10^® ® eV. The 
left panel shows the correlation function normalized to the phase space. The dotted 
line represents the expected distribution for an isotropic flux. The right panel is the 
corresponding statistical signihcance of the deviation from the background. 
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Figure 4: Signihcance of the self-correlation for published AGASA data with E > 10^^ ® 
eV relative to an isotropic hux at Earth. The left panel is the signihcance of the 
symmetrized correlation and the right panel shows the energy-ordered signihcance for the 
same events. Note that ID and 2D self-correlation plots, similar to those showed here, 
were hrst presented by the AGASA collaboration at the IGRG 2001 M. 


Assuming that this correlation is due to point sources, i.e. several of the cosmic rays 
originate from the same source, one can interpret the energy-ordered plot as follows. The 
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Figure 5: Sketch of two proton trajectories and the corresponding magnetic held di¬ 
rection. For the depicted conhguration of Ei > E 2 one gets hi < 62 independent of the 
latitude of the source. 

highest energy cosmic rays of the sample propagate with the least dehection. For negative 
A 6 this means the lower energy cosmic rays from the same source appear to come from a 
higher Galactic latitude (see Eq. (|)). This is only possible if they get dehected towards 
the Galactic south relative to the path of the higher energy cosmic rays. Depending on the 
sign of the charge this deflection determines the effective direction of the GMF transverse 
to the cosmic ray trajectory. However, only the effective direction integrated along the 
trajectory can be derived. For example, the GMF has to have a component parallel to 
the j/-axis for a positively charged particle from the source direction I k, 180°, as shown 
in Fig. 1^. The absolute value of the displacement of the correlation maximum from the 
center (A/ = A6 = 0) is given by the strength of the magnetic held and the energy 
difference between the cosmic rays comprising pairs with small AZ and A6. Of course, 
if the number of cosmic rays with positive and negative charge are approximately equal, 
one expects a symmetric signihcance pattern even for the energy-ordered correlation plot. 

We will discuss the AGASA data, including three of the observed multiplets, in more 
detail in Sec. ||. 
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4 Uniformly distributed sources 

The first scenario we explore is that of uniformly distributed sources. We inject protons of 
energy above 10^® eV from sources distributed uniformly on a sphere around the Galactic 
center and of radius Vgrc- In the simulation each cosmic ray is injected from a new source 
position. Astrophysically this would, depending on the value of rsrc, correspond to either 
sources in the Galactic halo, or a uniform and homogeneous distribution of extragalactic 
UHEGR sources in the absence of large-scale extragalactic magnetic helds. Here it is 
assumed that each source injects a very small flux, contributing not more than one particle 
to the cosmic ray spectrum seen at Earth. 

Even in the absence of the Galactic magnetic held one does not expect a completely 
isotropic hux at Earth. For isotropically emitting sources with locations symmetric about 
the Galactic center one does hnd an enhanced hux from the directions I = 90° and 270°. 
The enhancement is caused by the ohset of the position of the Earth from the center of the 
sphere where the sources are distributed (i.e. from the Galactic center). The enhancement 
occurs because for any longitude but I = 0 ° and I = 180°, the 1 /r^ behavior of the hux 
does not fully compensate the increase in the area (as seen from Earth) where the sources 
are located (see, for example, the discussion in p^). The enhancement factor is given by 


- 1/2 



( 10 ) 


Here i ?0 = 8.5 kpc is the Galactocentric distance of the position of the Solar System. For 
Galactocentric source distances of rg^c = 15, 20 and 40 kpc the maximum enhancement 
factor is ~ 1.214,1.105 and 1.023 respectively. 

Neglecting for the moment this small hux enhancement, one expects an isotropic hux 
at Earth even in the presence of the Galactic magnetic held. This is shown in the left- 
hand panel of Fig. and it is an ehect of Liouville’s theorem. The magnetic held does not 
change the magnitude of the cosmic ray momentum. Hence along each particle trajectory 
the phase space is preserved and the hux constant. If there exists a trajectory for any 
direction at Earth which can be followed back to the outer sphere of sources, the hux at 
Earth has to be the same as it is just outside the Galactic magnetosphere. If the hux is 
isotropic outside the GMF then it is also isotropic at Earth. 

Nonetheless it is interesting to see that the mapping of the sources to the arrival 
directions at Earth is realized in a non-trivial way. The sources that actually give rise to 
the hux at Earth are clustered in certain directions as shown in the right-hand panel of 
Fig. despite the fact source positions were drawn from an uniform distribution. The 
mapping reveals that sources in the Galactic plane are less efficient in producing GRs 
arriving at Earth. The reason for this can be understood with the aid of the right panel 
of Fig. 1^, which shows the projection onto the xz plane of three GR trajectories arriving 
at the detector from the same source in the Galactic plane. GRs have to travel along 
different magnetic arms in the Galactic disk in order to arrive at the detector. Each time 
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Figure 6: Left panel; equal area sky map of the arrival directions (in Galactic coordinates) 
of the cosmic rays with E > 10^®'^ eV. The sources from which protons are injected are 
uniformly distributed in a sphere centered around the Galactic center and radius r^rc = 40 
kpc. Right panel: sky map of of the position of the sources actually contributing to cosmic 
rays that arrive at Earth. 

an arm is crossed the held reverses sign dehecting the GR in opposite directions. Earth 
can only be reached by tracks that cross the Galactic plane at roughly the same location 
in between the source and the detector. This reduces the number of potentially detectable 
tracks compared to other source locations. 

Fig. I also shows a dehcit of sources in the direction of the Galactic center. The 
explanation of this effect is the large strength of the dipole held near the center of the 
Galaxy which strongly dehects the particles. 

These are two important features of the GMF model described above which will be 
rehected in all the scenarios we explore in this paper. In a model with diherent Galactic 
halo held, or diherent ratio between the strengths of the BSS and dipole helds, these 
ehects would change. 

Fig.0 shows another interesting feature of the GMF which is a tendency to dehect 
the GRs in the direction from the northern Galactic hemisphere towards the Galactic 
South. This ehect is caused by dehection in the BSS held, which is considerably stronger 
than the large-scale dipole held of the halo. The North to South dehection is caused by 
the direction of the local BSS held. This dehection is similar, although of much smaller 
magnitude, to the dehection in the Parker type halo held considered by Biermann et al. 

m- 

As is expected from Fig. P, the one-dimensional correlation function shown in Fig. is 
fully consistent with the expectations from an isotropic hux. Assuming the self-correlation 
observed by AGASA is not a statistical huctuation, a uniform and isotropic distribution 
of low-hux sources in the Galactic halo is not supported by current observations. 
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Deflection E>10 


19.6 



Figure 7: Sky map of the deflections of the cosmic rays that arrive at the detector. The 
arrows join the position of the source and the arrival direction of the cosmic rays coming 
from it. Only cosmic rays with E > 10^® ® eV are shown. The sources are located at a 
sphere of 20 kpc radius. 



Aa (deg.) Aa (deg.) 


Figure 8: Self correlation of arrival directions of MC generated events (uniform distri¬ 
bution of sources rg^c = 40 kpc) with E > 10^® ® eV. Left panel: correlation function 
normalized by the phase space. The self correlation for an isotropic background is also 
shown. Right panel: signihcance of the correlation (i.e. [signal-background]/width). 
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5 Supergalactic plane 


It is a long-standing question whether there is a correlation of the highest energy cosmic 
rays with the local matter enhancement, the so called supergalactic plane (SGP). An 

H 


indication for a large scale correlation was hrst reported in 


41 


The signihcance of 
and was replaced 


such a correlation decreased with the increasing world statistics 
by evidence for small scale clustering. A recent analysis of the world data on UHECRs 
concludes that the chance probability of the observed clustering within ±10° off the SGP 
is less than 1% [p. Furthermore the triplet and one of the two doublets in the AGASA 
data passing the outer Galaxy cut are located at less than 1° off the SGP. 
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Arrival directions E>10 ‘ eV 



Figure 9: Fqual area sky map (in Galactic coordinates) of the arrival directions of the 
cosmic rays with energy above 10^®'^ eV for Vgrc = 20 kpc. The sources are concentrated 
in the supergalactic plane (SGP). The solid line in the center marks the position of the 
SGP, whereas the left and right line around it bound the location where the sources are 
uniformly distributed (±20° in supergalactic latitude). 

To simulate sources in the SGP we sample source locations uniformly distributed in the 
direction of the SGP with a maximum offset of 20° in supergalactic latitude. Fig. ^ shows 
the corresponding arrival direction of protons emitted from sources at a Galactocentric 
distance of 20 kpc and E > 10^® '^ eV. Again we consider here the low-flux limit: each 
source emits only one particle. The total source efficiency is a convolution of the area 
within 20° off the supergalactic plane and the geometric effect of the right-hand panel of 
Fig. |. 
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Figure 10: Self-correlation of arrival directions of MC generated events with E > 
eV. The sources are distributed within ±20° off the SGP at a Galactocentric distance 
Tgrc = 20 kpc. Left panel; phase-space normalized correlation function. The self cor¬ 
relation for isotropic background events is also shown. Right panel: signihcance of the 
correlation in units of a. 

The mapping of the supergalactic plane to Galactic coordinates leads to a strong 
correlation in Galactic longitude. However the ID correlation in Fig. 0 is too broad as 
compared to data. The signihcance of the correlation is larger than 4a in an angular region 
of width ~ 30° which approximately corresponds to the band around the SGP where the 
sources are distributed. In addition, the features of the GMF already emphasized in the 
uniform source scenario are also seen. Gosmic rays are preferentially dehected to the south 
and certain source positions have a higher efficiency in producing observable cosmic rays. 
Both effects lead also to a very broad maximum in the two-dimensional correlation with 
a signihcance of up to 8a. 

It is interesting to note that the diherence in the real direction to the source and the 
arrival direction is large in latitude but not in longitude. Gonsequently the plane structure 
of the SGP in Galactic coordinates is approximately preserved by the GMF and the hux 
observed at Earth does not appear to be isotropic. 

We conclude that the SGP is disfavoured as the source of the events seen by AGASA, 
assuming a weak extragalactic magnetic held. The expected ID correlation could come 
closer to the observed correlation if the distribution of sources was restricted to a narrower 
band, say ±10° oh the SGP. However, the resulting arrival direction distribution would 
then be highly anisotropic, following the direction of the supergalactic plane. It can be 
seen in Fig. that the region around Galactic latitude 60°, where AGASA has seen several 
events, is largely unpopulated in this scenario even for sources within 20° oh the SGP. 
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6 Single point source 


The apparently isotropic arrival distribution of UHECRs strongly disfavours any scenario 
of a nearby, dominant single point source, unless a strong coherent large-scale magnetic 
held is assumed . 


The purpose of considering here single point sources is to study the correlation ex¬ 
pected from different source locations and source distances. The results are directly 
applicable to the interpretation of multiplets seen in data, under the assumption that 
they originate from point sources. 


In Fig. we show the signihcance plots of the symmetrized and energy-ordered self¬ 
correlation for single sources at four different directions. As discussed before, although the 
protons were only injected in a 30° cone around the line-of-sight to the Solar system the 
results apply to isotropically emitting point sources. Again the background is assumed 
to be an isotropic distribution, which leads to a large signihcance of the self-correlation. 
Only a number of representative examples of source locations in the AGASA outer Galaxy 
window are shown. 

The shape of the 2D correlations is strongly dependent on the position of the source, 
consistent with the dehections shown in Fig. ^ Due to the strong spiral component of 
the magnetic held parallel to the Galactic plane, the 2D correlation is more stretched in 
Ab than Al. Gosmic rays coming from sources with longitudes close to 90° get dehected 
toward larger (smaller) longitudes in the northern (southern) Galactic hemisphere. For 
sources in the vicinity of I = 180° there is almost no dehection in longitude. 

The dipole held would always dehect positively charged cosmic rays to larger Galactic 
longitudes. For the arrival direction this means that the cosmic rays appear to come from 
smaller Galactic longitudes (see Fig. ||). However, the small held strength of the dipole 
component makes it rather unimportant for the correlations considered here. Of course, 
the situation is diherent for cosmic ray trajectories coming close to the Galactic center, 
which are excluded due to the “outer Galaxy” cut. In particular for sources at high 
latitudes the non-trivial variation of the magnitude of the spiral held leads to a drift-like 
dehection parallel to the Galactic plane. 

The energy-ordered correlation plots show that, for positively charged multiplets of 
cosmic rays, the particle with lower energy is dehected more to the Galactic north and 
hence appears to come from further south (e.g. higher Galactic latitude in the southern 
hemisphere). The reason is the GMF component, parallel to the Galactic plane and 
pointing from large to small longitudes (see Fig. |I|). Since there is good agreement between 
diherent measurements on the direction of the GMF in the vicinity of the Solar system 


311, this prediction can be used to derive the charge sign of UHFGRs if point sources are 


identihed. 

A number of UHEGR models is based on point sources emitting cosmic rays only for a 
short period of time (e.g. Z-bursts, decay of super-massive particles). In Fig. we show 
the arrival time delay of cosmic rays coming from sources at a Galactocentric distance 
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of 20 kpc. The arrival time delay is defined as the time difference between the arrival of 
the cosmic ray and a light signal emitted from the source at the same time. For a given 
source, a strong correlation between energy and time delay is observed. Therefore one 
can estimate the expected time difference of cosmic rays emitted simultaneously from a 
point source simply by using this energy-time delay relation. Of course, such a simple 
energy-time delay relation is expected to hold only in the limit of small deflections. In 
this limit there exists essentially only one trajectory, up to small scale differences due to 
the local magnetic held huctuations, which connects the source with Earth for a given 
energy. 

The width of the distributions shows how strong the correlation is. The correlation 
will be wider if the fraction of the held strength coming from the turbulent component is 
increased. The simulations conhrm the expectation that high-energy cosmic rays arrive 
earlier than low-energy ones. However exceptions are possible. The width of the correla¬ 
tion allows for the inverted arrival time behaviour as long as the energies of the cosmic 
rays are diherent by less than 30%. 

Similarly, the angular diherence between the arrival direction and the line-of-sight 
to the source. Act, is directly related to the energy. With the 2D correlation plots and 
the knowledge of this angle one can estimate the source distance for a given multiplet 
and Galactic held conhguration. Fig. shows the expected angle Act for a number of 
source directions. The correlations depend only weakly on the longitude for hxed-latitude 
sources. 

Also clearly seen is the selection effect of the magnetic held. At certain energies 
almost no cosmic rays arrive at Earth form a given source. This can be studied in detail 
by comparing the injection spectrum with that of observed cosmic rays, see Fig. |T^. The 
low-energy hux from sources in the Galactic plane is strongly suppressed. On the other 
hand, the hux from sources at high Galactic latitudes does not suher energy-dependent 
modulation. 

Another way of characterizing the selection ehect due to the GMF is the source ef- 
hciency, dehned here as the ratio of the hux of detected GRs to the expected hux for 
vanishing GMF. Fig. shows the efficiency of a point source for diherent positions. The 
magnetic held has a focusing or defocusing ehect on the hux of cosmic rays emitted from 
a source depending on its position. Sources in the Galactic plane are less efficient than 
sources away from it as explained before. It is also interesting to note the slight asym¬ 
metry in source efficiency between the northern and southern Galactic hemispheres. For 
a hxed source position, the efficiency is a non-monotonous function of the energy. At 
sufficiently high energy the dehection in the magnetic held is small and the efficiency 
approaches 100%. 

So far we have considered only sources at a Galactocentric distance of 20 kpc and a 
Galactic magnetic held with a local strength of 1.5 fiG. We have checked with simulations 
for Tscr = 15, 20 and 40 kpc that the arrival time delay as well as the dehection angle 
depend only weakly on the distance to the source. However, as expected, the correlation 
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of both quantities with the arrival energy becomes weaker for increasing source distance. 
In contrast to the distance-dependence the amplitude of the GMF is very important. 
Fig. ^ shows the arrival time delay and deflection angle for a strength of 1.5 and 3 /xG of 
the regular component of the magnetic held in the vicinity of the Solar system. Doubling 
the held strength leads to an approximately twice as large dehection angle. The time 
delay increases by more than a factor of 2. 

7 Many point sources 

In the following we consider the case of many point sources. The hux of each of these point 
sources is large, i.e. each of them might contribute up to several cosmic rays observed 
at Earth. In astrophysical terms this scenario realizes some features expected in models 
for UHEGR production in the Galactic halo, such as Z-bursts, decay or annihilation of 
super-heavy particles, or acceleration by rapidly spinning magnetars. 

In the simulation the number of cosmic rays arriving at Earth from a particular source 
depends on the source location (i.e. its efficiency) and on the total number of cosmic 
rays injected per source. The latter parameter is one of the major unknowns in the 
simulation. In the following we assume that each source emits the same number of cosmic 
rays isotropically with a differential energy spectrum dN/dE oc i.e. all sources have 

the same cosmic ray luminosity. Thus, by construction, one has fewer cosmic rays injected 
from the same source at high energy than at low energy. This leads to a maximum of 
the self-correlation at some intermediate energy. At low energy (~ 10^® eV) the Galactic 
magnetic held will destroy the correlation and at high energy (~ 10^° eV) the number of 
cosmic rays from the same source is small. 

In the limit of many sources, the multiple point source scenario results in an almost 
perfectly isotropic arrival distribution on large angular scales, as already discussed in the 
uniform source distribution scenario. At the same time it can exhibit a signihcant self¬ 
correlation at small scales. Due to the large-scale isotropy we do not expect any signihcant 
diherences in the arrival hux for positively and negatively charged particles. Although 
diherent sources contribute to the hux at Earth, the observed hux will appear to be nearly 
isotropic. In particular there will be no asymmetry between the northern and southern 
Galactic hemispheres. 

In practice we inject 20,000 cosmic rays per source within a 30°-cone around the line- 
of-sight. The source positions are again sampled from a uniform distribution on spheres 
around the Galactic center. We do not restrict the position of the source, although for 
the analysis we only select events that come from the outer Galaxy as previously dehned. 
Simulations were done for the energy cutohs E > 10^® and eV. All cosmic rays 

arriving at the detector were recorded. The signihcance of the self-correlation in these 
data sets is then reduced by randomly diluting the samples of cosmic rays accepted for 
the analysis. For example, a set is diluted by a factor of 2 by accepting cosmic rays with 
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a probability of 0.5. 

In the following we consider two examples of typical correlations obtained in the 
multiple point source setup. Since the number of doublets depends on the total number of 
detected cosmic rays, we dilute our simulated data sets until we obtain a two-dimensional 
self-correlation of about 5-7 a at Al = Ab = 0. 

The hrst example shows the self-correlation plots for 200 cosmic rays arriving at Earth 
with energy E > 10^®'^ eV. The sources were located at a sphere with 20 kpc radius. The 
ID angular self-correlation in Fig. 0 exhibits a structure which is qualitatively similar 
to that seen in AGASA data: A large peak at small angular separation, followed by a 
tail at large angles. The signihcance plot reveals considerable statistical fluctuations due 
to the limited number of 200 cosmic rays used for this calculation. The observed cosmic 
rays were produced by 139 different sources, i.e. about 1.4 particles per source reached 
the detector. 

The 2D correlation calculated from the same events is shown in Fig. The structures 
of the correlation ellipse are qualitatively very similar to the one observed in the AGASA 
data (see Fig. ||). This could be considered as a strong support for our magnetic held 
model and point sources at 20 kpc Galactocentric distance. 

However, it should be emphasized that with only 200 cosmic rays the slope as well as 
the lobe-structures of the correlation ellipse are not statistical signihcant. To show this we 
calculate the self-correlation for 400 events, including the 200 previously analysed events, 
in Fig. 1^. The 400 cosmic rays at Earth were produced by 280 sources which again 
corresponds to about 1.4 detected cosmic rays per source. Although the self-correlation 
seems to be somewhat weaker and the slope of the correlation ellipse has changed, the 
signihcance of the excess at small Al and Ah is still obvious. 

Furthermore we want to emphasize that the energy-ordered correlation is always asym¬ 
metric towards positive Ah. This means that low-energy particles appear to arrive from 
directions with smaller b than high-energy ones. Such a result is expected from the held 
direction of the magnetic arm next to our Solar system. However, there is no such indi¬ 
cation in the AGASA data. In case of a mixture of protons and antiprotons as UHEGRs, 
as expected from decays of super-massive particles or neutrino annihilation, the energy- 
ordered correlation function is nearly symmetric. Our simulations show that this expected 
symmetry is only realized for samples with very high statistics. 

Finally we show in Fig. ^ the arrival distribution of the 200 cosmic rays used previ¬ 
ously in Figs. 0 and 1^. The location of the sources is also shown as the limited number 
of sources might contribute to some anisotropy. A visual inspection of Fig. reveals 
no signihcant large-scale anisotropy. We have also performed simulations restricted to 25 
cosmic rays arriving from the outer Galaxy, the same number as the AGASA data. The 
arrival distributions appeared to be compatible with that expected for an isotropic hux 
on large angular scales. 
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Point source at I = 112.5°, b = 0°: 



Point source at I = 112.5°, b = 45°: 



Point source at I = 112.5°, b = —45°: 



Point source at I = 157.5°, b = 45°: 
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Figure 11: Symmetrized (left) and energy-ordered (right) significance of self-correlation 
for E > 10^®'® eV. The plots were made from 200 cosmic rays arriving at Earth from 
sources with a Galactocentric distance of Tgrc = 20 kpc. 
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Figure 12: Arrival time delay of cosmic rays propagating from sources to Earth. The 
sources are again at a Galactocentric distance of Vgrc = 20 kpc. Each simulated cosmic 
ray trajectory, arriving at Earth from the “outer Galaxy” direction, is represented by one 
point in this plot. Gosmic rays were injected with a dN/dE ~ E~^'^ spectrum to obtain 
good statistics at high energy. 



Figure 13: Angle between arrival and source direction Aa for cosmic rays propagating 
from sources at a distance of 20 kpc from the Galactic center. 
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Figure 14: Injection and arrival spectra of cosmic rays for two source positions. The 
injection spectrum is proportional to The normalization is arbitrary. 




latitude b (deg.) latitude b (deg.) 


Figure 15: Efficiency (in %) of a single source for different source locations in Galactic 
coordinates and r^j.^ = 20 kpc. In the left panel the efficiency is calculated for cosmic rays 
with E > 10^®'*’ eV and in the right panel for cosmic rays with E > 10^®'^. The differential 
injection spectrum is proportional to 
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Figure 16: Arrival time delay and deflection angle for two different magnetic field 

strengths at the position of the Solar system. 
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Figure 17: One-dimensional self-correlation for events with E > 10^®'^ eV in a many 
point source scenario. The correlation is calculated from a data set containing 200 cosmic 
rays observed at Earth from sources at 20 kpc from the Galactic center. 
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Figure 18: Symmetrized (left) and energy-ordered (right) significance of self-correlation 
for E > eV. The sources are located at a r^rc = 20 kpc sphere around the Galactic 
center. The 200 events used for the calculation are the same as in Fig. 0- 
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Figure 19: Symmetrized (left) and energy-ordered (right) signihcance of self-correlation 
for E > 10^®'"^ eV. The sources are located at a = 20 kpc sphere around the Galactic 
center. The correlation was calculated from 400 trajectories reaching Earth. 
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Figure 20; Source and arrival direction distribution for 200 simulated cosmic rays with 
energy E > 10^®'^ eV. The injection spectrum at the sources was dN/dE oc E~‘^''^. 
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8 Application to AGASA multiplets 


As an application of onr resnlts we now discnss whether the clnsters fonnd in the AGASA 
data above 10^® ® eV are compatible with the model of the GMF we have adopted and 
the hypothesis of single point sonrces for their origin. Fig. shows the fnll data set with 
2 donblets and 1 triplet from the “onter Galaxy” . 
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Fignre 21: Arrival direction distribntion of cosmic rays with E > eV as pnblished 

by AGASA |^. The three mnltiplets listed in Table ^ are marked by sqnares (see text). 


For simplicity we nse only protons in onr simnlations. The resnlts do not change 
mnch for antiprotons bnt wonld be considerably different for nnclei. The approximate 
Galactic coordinates of the potential sonrce position are estimated based on the energy 
of the constitnents of the mnltiplet and the expected angnlar deflection of their tracks in 
the GMF. Then we inject protons from these potential sonrces and compnte, for different 
distances to the source, the time delay At between the arrival at the detector of a proton 
and a light ray emitted in the same direction, and the angular difference Aa between 
the line of sight from the detector to the source and the proton’s arrival direction. Our 
results are summarized in Table in which we also include the energy, arrival directions 
and dates of detection of events belonging to clusters as measured by AGASA . At and 


Aa depend on the energy of the events, the position of the source and the injection angle 
from the source. Protons having the same energy encounter slightly different magnetic 
held conhgurations along their trajectories to the detector depending on the injection 
angle. This produces the spread in At and Aa shown in Table ^ At and Aa are 
expected to have a maximum uncertainty due to the size of the detector of the order of 
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0.8 (rsi.c/20 kpc) years and 0.5° respectively. 


Table 1: Energy, Galactic latitude (6) and longitude (/) and date of arrival of the events 
belonging to clusters from the outer Galaxy as detected by AGASA. Rows from 9 to 12 
are the minimum and maximum values of the time delay between the arrival of a proton 
and a photon at the detector for sources at Galactocentric distances Tstc = 12, 15, 20 and 
40 kpc. Act is the angular difference between the line of sight from the detector to the 
source and the proton’s arrival direction. It is also computed for the same Galactocentric 
distances to the source as At. r® is the distance from the potential source of the multiplet 
to the Earth, calculated as well for the four Galactocentric distances. 


Cluster 

Cl 

C2 

C4 

E (eV) 

2.13 X 10^0 

5.07 X 10^® 

7.76 X 1019 

5.50 X 1019 

5.35 X 1019 

5.47 X 1019 

4.89 X 1019 

b 

-41.4° 

-42.5° 

55.1° 

56.6° 

56.2° 

-11.2° 

-10.8° 

1 

130.5° 

130.2° 

145.5° 

143.2° 

147.5° 

170.4° 

171.1° 

yr/mo/da 

93/12/03 

95/10/29 

95/01/26 

92/08/01 

98/04/04 

86/01/05 

95/11/15 

r® (kpc) 

5.2 


5.6 


3.5 

c© (kpc) 

8 

.8 


9.2 


6.5 

c© (kpc) 

14.4 


14.8 


11.5 

4° (kpc) 

35.1 


35.6 


31.5 

Ati 2 (yr) 

0.6-0.7 

12.5-15.8 

5.7-8.0 

11.0- 15.4 

11.7-17.0 

13.6-17.8 

16.7-23.3 

Ati5 (yr) 

1.4-2.1 

24.4-33.4 

11.8-16.0 

22.1 - 31.8 

22.3 - 33.0 

30.1-51.6 

37.1 -67.3 

At 20 (yr) 

1.6-2.7 

29.9-38.1 

14.2-20.8 

30.2 - 48.3 

31.4-48.3 

64.2- 77.1 

72.7-93.0 

At4o (yr) 

1.2-2.2 

32.4- 47.0 

16.8-28.6 

34.3 - 55.7 

41.2-63.9 

79.7- 108.2 

93.8 - 135.7 

Aai 2 (deg) 

0.4- 1.6 

3.5-5.5 

2.5-3.7 

3.6-5.0 

3.9-5.3 

5.2-6.9 

5.2-7.2 

Aai 5 (deg) 

1.1-2.0 

5.6- 7.1 

3.8-4.7 

5.1 - 6.6 

5.1 - 6.9 

6.6-9.5 

7.5 - 10.3 

Aa 2 o (deg) 

1.6- 2.7 

6.2-8.0 

4.0-5.1 

5.9-7.8 

6.1 - 7.8 

7.3-9.7 

7.9-10.8 

Aa4o (deg) 

1.2- 1.9 

6.3-8.8 

4.5-6.1 

6.6-8.5 

6.9-9.0 

5.5-7.9 

5.5-8.5 


As already discussed in Sec. H, due to the spread in At, the lowest energy event in 
the cluster does not necessarily arrive after the highest energy event and vice versa, as 
it would be expected from the larger angular deflection in the GMF of the lower energy 
event. For this to happen the energy of the events has to be roughly within ~ 30% of 
each other as can be understood by looking at Fig. |T^. For the sake of illustration of 
this point, take for instance the G4 doublet in Table |l| and rg^c = 15 kpc. Within the 
model of the GMF we have here adopted, the low energy event in the doublet could arrive 
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~ 40 years after a photon emitted in the same direction whereas for the high energy 
event this same delay could be ~ 50 years. Then assuming they are emitted at the same 
time, the high energy event would arrive 10 years after the low energy one. This is of 
course not compatible with the difference in arrival time observed by AGASA, however it 
is still possible and illustrates that one should be careful when ruling out a simultaneous 
emission from a source on the basis of the arrival times at the detector. 

For a cluster to be explainable as coming from a bursting source within our GMF 
model, the predicted values of At must be such that the time difference between the arrival 
of the events in the cluster is of the order of that observed by AGASA. Their angular 
separation at the detector should also be smaller than AGASA’s angular resolution of 
2.5°. With these two criteria in mind, the doublet G1 cannot be produced by simultaneous 
emission from a source within the range of Galactocentric distances 12 - 40 kpc, since 
the time between the arrival of both GRs in the doublet (~ 1.9 yr) is always smaller than 
the time difference predicted by the model, even when the spread in At is taken into 
account. This means that either the GMF strength is smaller than the one we are using 
by roughly an order of magnitude, or the source emits continuously and then the GRs are 
not emitted at the same time. It is also possible that the source is at a Galactocentric 
distance smaller than 12 kpc. In the G4 doublet, the two GRs arrive ~ 9.8 years apart 
in time and they can be produced by a bursting source for all source distances we have 
explored. The angular difference between the event’s tracks is smaller than 2.5°. The 
highest energy member of the G3 triplet arrives earlier than the event next in energy, 
but this arrival time inversion cannot be explained when accounting for the spread in At 
unless Tgj-c < 15 kpc. This points to a continuous source for the origin of the triplet, or to 
a nearby short-lived source, assuming the GMF model is correct. 

Finally we show in Table § the same information as in Table Q for three hypothetical 
clusters of events extracted from the AGASA’s sample. They are constituted of events 
which are close to each other in the sky as can be seen in Fig.|^. Their angular separation 
is larger than 2.5°, being ~ 5.6° for the events in the D1 pair, ~ 5.3° for the events in 
D3 and from ~ 3.3° to ~ 7.4° for the three events in D2. However, the arrival directions 
follow the expected pattern of deflection due to the GMF (see Fig. |^. 

We explore the possibility of each of them being produced by a distinct point source 
on the basis of their angular separation. Due to the spread in Aa the two members 
of the D1 hypothetical doublet could qualify as a doublet. The GMF model we use 
predicts that the lower energy protons from a cluster should come from smaller values of 
b. The opposite is observed for the pair of events that constitute the D1 doublet, however 
the difference is still within the experimental uncertainty in the determination of arrival 
direction. The range of angular separation between the members of the potential triplet 
D2 would in principle allow that they belong to a triplet. However the sign of Ab is the 
opposite of what is predicted in the GMF model and, unlike doublet Dl, is much too 
large to be attributable to the experimental uncertainty. A possibility which is not yet 
excluded is that the low energy member of D2 is an antiproton which would explain why 


Table 2: Energy, Galactic latitude (6) and longitude (/) and date of arrival of the events 
detected by AGASA from the outer Galaxy that could belong to hypothetical clusters. 
Rows from 5 to 8 are the minimum and maximum values of the time delay between the 
arrival of a proton and a photon at the detector for sources at Galactocentric distances 
Tsrc = 12, 15, 20 and 40 kpc. Aa is the angular difference between the line of sight from the 
detector to the source and the proton’s arrival direction. It is also computed for the same 
Galactocentric distances to the source as At. re is the distance from the potential source 
of the multiplet to the Earth, calculated as well for the four Galactocentric distances. 


Cluster 

D1 

D2 

D3 

E (eV) 

6.81 X IQi® 

4.95 X 1019 

7.16 X IQi® 

6.11 X IQi® 

4.29 X IQi® 

9.10 X 10^9 

4.78 X 1019 

b 

-15.7° 

-14.0° 

-4.0° 

-7.8° 

-0.7° 

25.6° 

29.8° 

1 

93.3° 

98.8° 

149.8° 

152.4° 

150.3° 

108.8° 

105.1° 

yr/mo/da 

84/12/12 

99/09/25 

99/07/28 

98/10/27 

99/10/20 

91/11/29 

96/05/13 

r® (kpc) 

7.1 


3.8 


6.2 

ce (kpc) 

11.0 


6.9 


9.9 

r|0 (kpc) 

16.7 


12.0 


15.6 

4° (kpc) 

37.6 


32.2 


36.4 

<1 

18.5-26.0 

37.3 - 52.9 

4.8- 7.5 

6.8-10.6 

14.8-21.5 

5.3-8.0 

16.4-29.8 

Ati5 (yr) 

15.6-28.5 

30.0-61.3 

13.7-23.2 

17.8-33.1 

49.0-86.2 

7.9-14.0 

29.7- 54.2 

At 20 (yr) 

36.1 -46.3 

62.7- 93.8 

12.8-28.4 

16.1-41.3 

37.6-93.0 

7.9-15.1 

31.5-57.6 

At40 (yr) 

77.4-254.1 

148.4-441.3 

25.1 -39.9 

35.3- 52.6 

63.0 - 105.0 

11.8-20.3 

37.9- 79.7 

Aai 2 (deg) 

3.6-4.6 

4.9-6.8 

3.0-4.4 

3.7- 5.0 

6.0- 7.1 

1.6-2.9 

3.4-5.4 

Aai5 (deg) 

4.1-4.9 

5.8-7.7 

5.1-6.9 

6.2-8.4 

9.5-12.7 

3.1 - 3.8 

5.8-6.9 

Aa 2 o (deg) 

3.5-5.1 

4.8-6.5 

4.0-6.3 

5.0-8.0 

8.0-12.3 

2.6-3.7 

5.5 - 7.6 

Aa4o (deg) 

0.8-4.7 

1.0-4.5 

2.7-5.4 

3.5-6.5 

4.9-9.8 

2.6-4.0 

3.8 - 8.9 


it is arriving from a larger b than the high energy members of the triplet. However if this 
is the case, the sources have to be nearby otherwise the angular difference would be too 
large. The same comment applies to the constituents of the hypothetical doublet D3. The 
D1 doublet is hence not yet ruled out as coming from a single source within the GMF 
model used in this paper. Glusters D2 and D3 are less favored but still possible in some 
scenarios. Future accumulation of statistics by AGASA and other experiments might well 
hud actual clusters of events close to the directions of the three potential clusters we have 
just discussed. 
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9 Conclusions 


We have performed a detailed calculation of the UHECR flux expected in several simple, 
limiting source scenarios. We have implemented an approach different from that used by 
other authors propagating the particles from the sources to the detector. This allows not 
only the investigation of the particle’s deflection in the GMF but also the detailed study of 
various source scenarios. We hnd the self-correlation analyses in one and two dimensions 
to be powerful tools for identihcation of the sources of UHECR and their distribution. 

The predictive power of our results is mainly limited by the uncertainty in the GMF 
model we use in terms of structure, held strength and spatial extent. Nevertheless we 
have obtained a number of interesting results which do not depend very much on the 
parameters of the particular GMF model. 

• For an isotropic extragalactic hux one expects an isotropic hux at Earth as long 
as cosmic rays are not trapped by the GMF. The observed CRs would stem from 
sources being distributed non-uniformly. 

• The SGP plane is disfavoured as main source of UHECR as long as the extragalactic 
magnetic held is weak. 

• Point sources produce two-dimensional correlation ellipses, the slope of which de¬ 
pends on the position of the source and not on its distance. 

• The observed cosmic ray hux of point sources is energy-dependent modulated. For 
example, cosmic rays with 10^® eV are suppressed from sources in the Galactic plane. 

• The size of the correlation in Al — Ab depends strongly on the strength of the GMF 
and only weakly on the distance to the sources. 

• The energy-ordered correlation analysis reveals the general held direction/polarity 
in the vicinity of the Solar system. It can be used to determine the charge sign of 
the cosmic rays and also helps to estimate the importance of statistical huctuations 
in data. 

• The width of the time delay versus energy correlation allows for the inverted arrival 
time behaviour, i.e. the higher energy cosmic ray of a pair arrives earlier than the 
lower energy one. This may happen if the energies of the cosmic rays are within 
~ 30% of each other. 

• It seems conceivable that a scenario of multiple, uniformly distributed sources can 
produce correlation patterns as seen in AGASA data. 

• The slope and the lobes of the 2D correlation ellipses are not statistically signihcant 
for less than ~ 200 events from the outer Galaxy in the multiple sources scenario. 
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• The energy-ordered 2D correlation distribution of the AGASA data does not agree 
with the expectations for positively charged particles and the BSS GMF model, 
however the small statistics of the data does not allow us to draw any hrm conclu¬ 
sions. 

• The analysis of the AGASA multiplets within a model of point sources shows that 
only one doublet is excluded from being a pair of cosmic rays simultaneously emitted 
from a source at a distance between 15 - 40 kpc. All other multiplets are within the 
uncertainties compatible with coming from short-lived point sources. 

• Further three possible multiplets were selected on the basis of the expected self¬ 
correlation patterns. 

In forthcoming work we plan to derive limits on UHEGR models such as Z-bursts and 
decay of super-heavy particles. Supposing the self-correlation observed in AGASA data 
is not a statistical fluctuation, we will extend our work to set limits on the minimum 
number of nucleon-antinucleon pairs needed per point source to produce doublets and 
triplets as well as limits on the total number of point sources required to sustain the 
observed UHEGR flux. 
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